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We discuss the problem of an effective descriptions of the phase transition phenomena 
' in the pure gluodynamics in SU(2) symmetric QCD. We choose the method of calculation 

following the conjecture that the infrared sector of the theory possesses the same confinement 
characteristic as the full theory. We show, that analytic descriptions of this phenomena 
f*) • is beyond the Gaussian method of evaluations of functional integrals. We propose a non- 

perturbative evaluation of functional integral, meanwhile for two dimensional Wiener integral 
for x theory. 
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INTRODUCTION 



The evaluations of the measurable quantities in quantum field theory (QFT) are well - defined 
theoretically, but in real calculations appear many obstacles connected with the non-ability to 
realize the calculations corresponding to the theoretical definitions. The calculations in QFT are 
usually performed analytically using perturbative calculus, numerically by lattice gauge theory 
(LGT), or by creating models based on the results of previous methods of calculations. The 
perturbative methods give analytical results, but various problems in QFT are non-perturbative 
phenomena, therefore they are applicable for the narrow band of problems only. LGT, the best 
method of evaluations in this time, suffers from finiteness of the lattice, problems of definitions of 
the mathematical objects on the lattice and the continuum limit problems. We adopt the method 
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of analytical calculations in the continuum to study of the transitions in dense matter. After short 
introduction into the problem in the next section we propose in the section 3 the non-perturbative 
solution of the problem of the evaluation of the functional integral with fourth order term in the 
action. We found the result in the form of an asymptotical series, and a method of solution of the 
problem is described in papers [lj and Qj. 



II. EFFECTIVE DESCRIPTION OF THE CONFINEMENT IN CONTINUUM 

The Euclidean finite temperature theory is defined by the functional integral for partition func- 
tion: 

Z((3) = N J [DAD ex P("^) 
where the Euclidean action is defined by: 

Se = \[ dT j d?>x:F ^ T % 
The Euclidean color field strength is defined by: 

•F \iv ®iJfAv &vA^ gs A^A U 

We require the periodicity of color potentials in the direction of the imaginary time variable: 

A%{T + /3,x)=AftT,x) . 

The theory with corresponding partition function looks like a partition for the finite temperature 
theory in statistical physics. This "isomorphism" allows us to investigate the QCD problem by 
methods of statistical physics. 

We have concerned on the phase transitions to the quark gluon plasma at high densities of 
the matter, the one of the challenging forecast of the QCD. We investigate some problems of the 
study of this phase transition, usually known as confinement /deconfinement phase transition, in 
model independent method described bellow. Such description of the confinement problem rely 
on the conjecture that the infrared sector of the finite temperature theory possesses the same 
phase structure as the complete theory. This is the assertion of Appelquist - Carrazone decoupling 
theorem [3] , where it was proven that integrating off massive modes of the theory one obtains terms 
in the potential controlling dynamics (in this case the infrared variables) of the theory. By help 
of this theorem we reduce the infinite number of the variables of the original theory to the finite 
number in the infrared sector effectively describing the phenomena studied. 
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To evaluate the integration over the massive modes we proceed as follows: 

1) We separate the field variable to the sum of the "classical" part and "quantum fluctuation": 

Al(t,x)=Al(t,x) + al(t,x). 

We use the static gauge for evaluations of the functional integrals over massive modes, The saddle 
points of the action in this gauge are fixed as: 

Af = 0, i,a= 1, 2, 3, A\ = A% = 0, A% = const. 

We choose finite ^4q in the color space for the simplicity. 

2) The periodicity of the color potentials in the imaginary time variable offers a chance to use 
Fourier transformation of the potentials in this variable. This operation automatically separates the 
Fourier transforms of the potentials to the infrared, zero frequency modes and massive modes, where 
nonzero fourier frequency modes (= Matsubara frequencies) appear as mass terms. Integrating over 
this "massive modes" by Gaussian method we obtain the contribution to the effective potential 
(in the scope of this article named V e ff = V e ff(AQ, 6q), 6°, a = i = 1,2,3) and contributions 
of the Aq dependent mass terms for the perpendicular quantum fluctuation degrees of freedom 
(bl,bj,i = 1,2,3) in the infrared sector of the theory. This process is known as dimensional 
reduction, described by Appelquist and Pisarski [4], offers possible calculation scheme for finite 
temperature field theories. The zero modes action in quantum fluctuations fields read: 

S zero = P f ^{Ig^G^ + |Cff^S) 2 [(6j) 2 + (6?) 2 ] + |(^) a + 6?)} (1) 



with field strength 



Gfj = dib] - djbl - ge abc b\b] 



We have found the two different effective systems, depending on the value of the static color 
potential Aq. 

1. Let A% = 0. In this case all quadratic terms in zero frequency field variables in Eq.([T]) 
disappears. The effective infrared system depend on variables a,i = 1,2,3. Color SU(2) 
symmetry of the system is maintained. This system correspond to interaction of the effective 
Higgs field, represented by chromoelectric potential 6q, with chromomagnetic fields represented by 
chromomagnetic potentials W. 

2. Let Aq 3> 0, at least as lowest Matsubara frequency. In this case the modes b\,bf in Eq. 
(PQ) must be treated as massive, and we integrate off these modes from the zero mode system. 
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Then the infrared system depend on variables 6^,6q, i = 1,2,3 and the remanent color symmetry 
group is Z(2). It is the well-known abelian projection observed in the lattice calculations also. In 
earlier continuum theory studies Q], Q], we find in this sector the confinement/deconfinement 
transition with correct behavior of the potential between heavy quarks calculated as the corellators 
of Polyakov's loops. 

But what happens when Aq grows from zero to a finite value and infrared system undergo 
transition SU(2) — > Z(2)? This phenomenon we „a m e d as "abe— — , Q. For 
its studies the use the Gaussian method for evaluations of functional integrals is not sufficient. 
Preliminary investigations indicate that this transition could be non-continuous, and consequently 
the lattice formalism may be not powerful also. We conclude, that one needs a formalism for the 
evaluations of the functional integrals with a higher power terms in the potential. This formalism 
enables to eliminate the mass term, even in this limit the functional integral is well defined. 

As a first attempt we applied the method to evaluate (semi) analytically Wiener functional 
integral with x 4 term in the action. The short description of this method is reported in the next 
section. 

III. EVALUATION OF ip 4 WIENER FUNCTIONAL INTEGRAL 

The simplest non-gaussian functional integral is Wiener functional integral with x 4 term in the 
action. In Euclidean sector of the theory we have to evaluate the continuum Wiener functional 



integral: 




(2) 



In this case the action possesses the fourth order term: 
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The continuum Wiener functional integral is defined by a formal limit: 



Z = lim Z 



■N ■ 



Here the finite dimensional integral Zjy is defined by the time 



slicing method: 



Zn 




(3) 
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where A = (3/N, and a,b,c are the parameters of the model. The quantity Zn represents the 
unconditional propagation from (p = <po = to any <p = ipjq (Eq. ([3]) contains an integration over 
<Pn)- 



A. One dimensional integral 



An important task is to calculate the one dimensional integral 

+oo 

h = J dx exp{-(^x 4 + Bx 2 + Cx)} , 



where ReA > 0. The standard perturbative procedure rely on Taylor's decomposition of exp(— Ax^) 
term with consecutive replacements of the integration and summation order. The integrals can be 
calculated, but the sum is divergent. 

Instead we propose the power expansion in C: 



-C) n 



n=0 



dx x n exp{-(Ax 4 + Bx 2 )} . 



The integrals in the above relation can be expressed by the parabolic cylinder functions 
D_ v _i/ 2 {z). Then, the integral I\ can be reds: 



r(l/2) f P ( w c 2 B 



B 1 ml 

m=U 



2A 



(4) 



where we used the abbreviation: 

V_ m _ 1/2 (z) = z m+1 / 2 eT D_ m _ 1/2 {z) . 

It was shown, that sum in the Eq. (U]) is convergent and for finite values z this sum converges 
uniformly in £. 

Applying this idea of integration to the iV dimensional integral ([3]) we proved [1] the exact 
formula: 



■>N 



N 



H 2(1 + 6A 2 /cVi 



i=0 



with 



oo N 

e n 

fci,— ,fcjv-i=0 «=0 



'2 k r 



(2fci)! 



r(/^_! + fci + 1/2) ^ x'-^.i-fei-irt 
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Fig. 1: b dependence of the continuum function S(a, b, c, r) for fixed values a, c, r. The first three nontrivial 
terms of the asymptotic series j5]) were used. 



where k = k N = 0, p = (I + &A 2 /c) -1 , z = c(l + 6A 2 /c)/ v / 2aA 3 , = 1 - vl 2 /wi_i, w = 
1/2 + A6A 2 /c, A = 2(1+ ; ) 1 A ij /e) • We see that 

- p is independent of the coupling constant. 

- Only the argument z of parabolic cylinder function is coupling dependent. 

To evaluate Sjv, we must solve the problem how to sum up the product of two parabolic cylinder 
functions. The parabolic cylinder functions are related to the representation of the group of the 
upper triangular matrices, so we implicitly expect a simplification of their product. This problem 
is not solved completely yet. We adopt less complex method of summation, namely we use the 
asymptotic expansion of one of them, then, exchanging the order of summations we can sum over 

Sure,, the result is degraded to an asymptotic expansion only, bnt still we have an analytical 
solution of the problem. This procedure was widely discussed in detail in paper I [1|], here we 
remind the result: 



fj.=0 v ' 

The evaluation of the symbols { C 2/i (N) } Q is described in paper II [2]. We explicitly present the 
first nontrivial contribution in the continuum limit (p finite, A — > 0, AiV = (3): 

3 

{C 2 (a,b,c,T)} 2Q = ^-g [37rtanh 2 (7r) + tanh^r) -77-] (6) 

where 7 = ^J2b/c. A typical b dependence of the S(a,b,c,r) = limTv^oo Sn is expressed on the 
Fig.l. 
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The parameters a, c, r are fixed constant, and we took J = 3 in ©. For 6 < we see the 
singularities for \ jt\ = (n + 1/2)tt. The divergent behavior correspond to the powers of tan (|7t|). 

B. Gel'fand - Yaglom Equation 

n 

Gel'fand and Yaglom proved in [9 ; ] for the harmonic oscillator that nontrivial continuum limit of 
the finite dimensional integral approximation to the functional integral should be evaluated from 
N— dimensional integral results by a recurrent procedure. For harmonic oscillator Gel'fand and 
Yaglom derived for the unconditional measure integral Z(j3) the differential equation: 



d 2 2b 

^r 2 y{r) = - y(r), 



(7) 



where Z(0) = lim^oo Z N = y^)" 1 / 2 

Following the idea of Gel'fand and Yaglom we found for an-harmonic oscillator the generalized 
Gel'fand - Yaglom equation. We define the unconditional measure functional integral Z{(5) by 
relation: 

i y(P) 



Z(f3) 



lim Zn 

N^oo 



, HP) 



The function S(r) is given as the continuum limit of Eq.([5 



S(t) 



lim 5/v ■ 

N^oo 



The generalized Gel'fand - Yaglom equation read: 



d 2 

0^y(r) = y(r) 



(8) 



accompanied by initial conditions y(0) = ^(O) 2 , and 



T = 



lS(r) 2 



T = 



For S(t) one can use a perturbative expansion in coupling constant a and then solve Eq. ((§]). This 
procedure gives a non-perturbative approximation of the functional integral (J2J. For harmonic 
oscillator limit we have S(t) — » 1. 

A typical dependence of the functions 4 ( In S(t)) 2 on r for positive fixed values of parameters 
a,b,c is shown on the figure 2. 



IV. CONCLUSIONS 



We stressed some features of the confinement problem in the continuum theory, which can be 
solved by non-perturbative method of evaluations of the functional integral. At the same time, in 
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Fig. 2: t dependence of the continuum function 4[9 P ln(5(a, b, c, t))] 2 for fixed a, 6, c. The first three nontrivial 
terms of the asymptotic series ([5]) were used. 

the second part of the article we reported a new approach to the solution of such problem. We find 
for the functional integral of an an-harmonic oscillator the non-perturbative equation ([8]) . Solving 
this equation we in principle can find the analytical solution of the an-harmonic oscillator problems 
as energy levels, measurable quantities, etc, the studies are in progress. 

The skeptic view of such solution of the problem may be expressed by the comment that we 
reformulate the problem to the language of differential equations. It is true. But our opinion is 
that theory of the differential equations is more elaborated and flexible than approaches based on 
naive perturbative theory and it can give more reliable results than the perturbative theory. 
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